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Abstract. The parametric instability of upper hybrid wave decay into back scattered 
upper hybrid wave and lower hybrid wave is considered for conditions of inhomogencous 
plasma of spherical tokamaks. The possibility of absolute instability is demonstrated 
and the corresponding threshold is determined. It is shown that the threshold 
power increases with pump frequency and electron temperature. Threshold power 
is estimated for typical parameters of experiment in MAST tokamak. It is shown that 
in this case parametrical reflection arises, if probing power exceeds 90 W/cm 2 , which 
gives 30 kW for a beam of 10 cm radius. 



1. Introduction 

In recent years, considerable attention of the controlled fusion community has been 
paid to spherical tokamaks (ST). These are small aspect ratio devices with typically 
high plasma density and comparatively low magnetic field. This ST feature has strong 
effect on the electromagnetic wave propagation. In the microwave frequency region, 
characteristic surfaces, like the upper hybrid resonance and the cut-off are very close to 
the plasma edge. As a result, the electromagnetic (EM) waves are unable to penetrate 
into the plasma interior. The only way to overcome this difficulty is to use the 
linear conversion of the incident EM wave into the electron Bernstein wave (EBW) 
occurring in the upper hybrid resonance (UHR). The latter has no density limitations 
and can, in principle, carry the radio frequency power deep into the plasma. This 
mechanism of wave conversion has been successfully demonstrated to produce heating 
in over dense plasmas in the W7-AS stellarator [3]. The plasma heating experiment 
based on this scheme is in progress now in the MAST tokamak at Culham, UK. The 
wave propagation, conversion in the UHR and absorption is usually accompanied, in 
100 kW power level experiments, by nonlinear effects, in particular, by parametric decay 
instabilities (Versator [T], FT-1 [2|, W7-AS [3]). These instabilities lead to redistribution 
of incident power between plasma species and can cause anomalous reflection, especially 
when excited at the plasma edge. 
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The present paper is devoted to analysis of the decay instability thresholds and 
growth rates for specific conditions of low magnetic field typical for ST. The study is 
focused upon the decay of UH wave into another UH wave and intermediate frequency 
range wave satisfying the lower hybrid resonance condition, which was observed in the 
UHR heating experiments mentioned above. The influence of plasma inhomogeneity on 
its threshold is investigated for backscattering of the incident UH wave. Dependence 
of the decay instability threshold on the pump frequency, necessary for the heating 
experiment optimization, is studied. 

The paper is organized as follows. In section |21 we deduce equations, describing 
the decay of the incident high-frequency UH wave into UH wave and low-frequency LH 
wave: £jjh — ► £'uh + ^lh, and consider them in WKB approximation. In section El we 
calculate an absolute instability threshold, which corresponds to UH wave induced back 
scattering instability. Brief discussion follows in section |3J 



2. Equations for wave amplitudes 

We use slab plasma model, i.e. density and magnetic field gradients are assumed to 
be along x axe. A magnetic field direction is chosen along z axe. We consider one 
dimensional problem of pump wave parametric decay. UH pump wave is supposed to be 
excited by external antenna via tunnelling effect (X— >B scheme according to [I]), and 
assumed to propagate in x direction. We consider here the high pump frequency case, 
when the frequency is larger than doubled electron cyclotron frequency, corresponding 
to the magnetic field in UHR: uoq > 2uj ce . In this case the UH pump wave dispersion 
curve (see figure does not possess a turning point and transformation to Bernstein 
wave occurs without change of group velocity sign. 

By indices o, i, 2 we will mark frequency, wavenumber, complex amplitudes of the 
electric fields and potentials of the pump wave, parametrically reflected UH wave and 
LH wave correspondingly. 

2.1. Nonlinear current and equation for LH wave 

Poisson equation for LH waves can be represented in the following form 



divD LH = 

dx 



^M+£ 2 tM—)e : 



LH 



+ rj(u 2 )—: = 47T p LH 

dz 



Here E LH = -V(j) LH is an electric field of LH wave, which is assumed to be potential, 
e, n are the components of the dielectric tensor 

£{U) = 1 2 ~ JL — 2 2 2~' = 1 2 

which for LH wave frequency uj 2 ~ ^uj ce u C i take the form 
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where Vr e ,i corresponds to electron and ion thermal velocity Vr e ,i = (2T e j/m^j) 1 ' 2 . In 
particular, for LH wave (J2J) takes the form 



3 w 



2 tot V4 



rrii 



4W- ^ x + ?^ 



Thus, equation describing the excitation of LH wave, can be rewritten as 

divD LH = -£|(u;2)0 / ^-e(cj2)0LH-e / (^2)0L/f + ^(^2)0LH = 4vrp LH (3) 

Here and below ' denotes d/dx. A charge density p is associated with nonlinear current 
3lh by continuity equation 



dpLH dj LH 







dt dx 

To obtain nonlinear current Jlh we consider electron motion in the field of three potential 
waves 

e 



v. 



2m P 



UJ rP V,, 



-ji^exp i / ki{x')dx' — iuJit + c.c. > 

>— n 19V L J J ) 



i=0,l,2 



U Ce Vy 



(4) 



Here dot ' means d/dt. Electric field of three waves is taken in geometrical optics 
(or WKB) approximation. This approximation is not valid if the decay point x d , 
determining in the inhomogeneous plasma by the conditions 

k (x d ) = ki(x d ) + k 2 (x d ), u = uj 1 + uj 2 (5) 

is situated in the vicinity of LH wave turning point (see section |U for proper discussion 
of corresponding criteria). 

In deducing (JTJ) we assumed following criteria to be satisfied 



^2^Te 



«1, 



k 2 z Vl 

UJn 



<1, 



&2 — &ce 



kM 



zVTe 



»1, 



uj 2 — 2u c 



z VTe 



> 1 



First criterion, which characterizes /^-approximation, allows us to get nonlinear 
component of a solution of (J3J) in the form 



VLH 



4m 2 



k; 



e ^=0,1,2^ ~< l(Wi + Wi) 2 -W, 



UJ; + UJ 



J 



EiEj exp 



i / (ki + kj)dx' — i(u>i + ujj)t 



UJ;, — OJj 



(Ui -LUj) 2 -Ldl 



-EiE*exp 



i / (fcj — kj)dx' — i(u>i — ujj)t 



c.c. 



Averaging vlh we neglect high-frequency terms. That yields 



(vlh) 



OJl - U! 



Am 2 (uj -u s ) 2 - u 2 ce \uo 2 - 



k n 



uJi 



E Q E{ 



x exp 



(ko — ki)dx' — i(oj — u>i)t 



c.c. 
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Taking into account that j LH = —en (v LH ) and passing to the complex amplitudes of 
the potential <Plh — (02 + ^2) /2 one obtains an equation for LH wave 



- 4(w 2 )C - - e'{u) 2 )<l> 2 + k 2 zV {^ 2 

,2 



e oj, 



pe 



fc fci(fco - h) 



2m e uj 2 e u iV 1; V — u 2 „ uj 2 — to 2 c 



fa 



fa 



x exp 



/X 
(fa — ki)dx' — i(cv — u>i)t 



(6) 



2.2. Nonlinear current and equation for UH wave 



For UH waves we have [3] 



D 



UH 



dx 



e(u ) + £ T (u ] 



dx 2 



E U X H + ig{u Q )E, 



UH 

y 



d2 E UH 



dx 2 



ig(u )-$-E : 



2 

rpUH 
2*x 



(7) 



Here corresponding components of dielectric tensor take the following form for the 
frequency of UH wave uj\ ~ u UH = ui 2 + to 2 e 



UJ 



pe 



u\ - u 2 e 



Parameter £t(wi) can be represented as 



VJ 2 

Te 



< £ 2 T M 



Considering potential UH wave E UH = —VtfiuH-, and substituting integrated second 
equation of (jZj) to the first, one obtains 



- l 2 T M<f>u H - e{^)<j>' UH - e'M^x + g 2 M^<f> 



UH 



UH 



(8) 



A charge density pun is associated with nonlinear current jun as 
dpuH , d(j UH ) 



dt 



dx 



0. 



Juh = -e8n W2 v u 



Here v Uo describes electron motion in the field of the pump wave 



lime u 2 - u 2 e 



iE exp 



k (x')dx' — iu t 



c.c. 



Density modulation Sn^ is caused by the electron motion in the field of LH wave 



1 ulM 



pe 



lire u 2 - u 2 e 



\E 2 exp 



i / k 2 (x')dx' — iuj 2 t 



c.c. 



Here we omitted a contribution of LH wave component along the magnetic field, which 
is smaller in factor of k 2 V^ e /u 2 <C 1. 

Averaging the nonlinear current, we leave the terms varying with frequency loq — uj 2 
only. That gives 

,2 



(Juh) 



e w 



pe 



u k 2 



16% m e uj 2 e (x! 2 



EqE 2 exp 



i / (fa — k 2 )dx' - i(uj - ui 2 )t 



c.c. 
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yielding (jSJ) for the complex amplitude of the potential 4>uh 



i>i + 0i)/2 in the form 



to 
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2.3. WKB-analysis of the equations obtained 

A dispersion relations, which can be obtained from equations 0, ©, when pun 
Plh = 0, take the following form pj: for UH waves 

ki 



e(o; ,i) = &r( w o,i) fc ,i 



k 2 



(10) 



where the transformation wavenumber is 

k 2 _ u} o g _ u ce /c 

and for LH waves 



c £ T (u ) £ t (u ) 
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Equations (fT0*|) , |TT|) allow us to obtain group velocities of the corresponding waves. We 
get 



<4l - ^le 
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(12) 



One can see from (|1() J) .(|11 J) . (j!2j) that in the probing frequency range under consideration 
ujq > 2uj ce , which is used at present for EBW heating in MAST, there is no change of 
the group velocity sign in the UHR point. The transformation point of LH wave, which 
is shifted from LH resonance position (where e(u 2 ) = 0), is the turning point of LH 
wave, and group velocity changes the sign there. Corresponding dispersion curves are 
represented in figure^ We consider the most interesting case of k Q > 0, k\ < 0, when the 
group velocity directions give rise to positive feedback loop, which can lead to absolute 
decay instability 0001111!. 

We consider equations ©, ® in WKB approximation, substituting 



a 0,l,2 



'0,1,2 



exp 



k 0jlt2 (x')dx' - iu 0t i, 2 t 



fco,l,2\/yo7l,2 

and neglecting corresponding small terms. In the vicinity of the decay point (0) we have 



v\ a 2 exp 
v 2 a\ exp 

,2 



U) 



pe 
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k 2 a* 



v 2 



to 2 



k 2 a 



4m e L)l - U 2 e y/-V ViV 2 



(13) 
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Figure 1. UH and LH waves dispersion curves in high frequency case (uio > 2ui ce ). 



3. Absolute instability threshold 



Absolute instability can arise, when decay conditions (0) allow two decay points x±, 2 to 
exist, and the group velocities directions provide positive feedback loop. In this case, 
according to [HI El; the absolute instability threshold is determined by the following 
conditions on the waves amplification coefficient 

\S 12 (x 1 )\\S 21 (x 2 )\ = l (14) 

where Sjk(xi) is the wave amplitude a^, which leaves the vicinity of decay point Xi, due 
to incidence onto this point of the wave aj of unit amplitude 



,5' 



12 



^\Z\/2 



r(i|z| + i)~ ~ 21 r(-i|z| 

1 2 v\v 2 and I is the length of the decay region 

^ -l 

£ 2 {xd) = { -j- [k (x) - fci(x) - k 2 (x)\ 



^n\Z\/2 



where Z 



dx 



x=x d 



The spectrum of the instabilities arising is determined by the condition on the 
phase, gained in the feedback loop 

[ifc (x') - h(x') - k 2 {x')] dx' + - = 2nN, N = 1,2,... (15) 



■n 



To calculate the decay instability threshold we act in accordance with following 
procedure. We calculate the terms, involved in (fH^ using (|TH|) . and then substitute 
them to (|T4*j) . obtaining an equation for threshold power. 

At first we calculate the decay point coordinates x± >2 . It should be noted that the 
UHR position xuni^i) of the parametrically reflected UH wave is shifted in respect to 
the UHR position of the pump wave xuh{^o)- This shift can be estimated as 

2L [x uh (uj )]uj uj 2 



xuh{oJi) - x UH {u ) 



L(x) 



1 dn(x) 
n dx 



^pe 

1 



pe 

2tul dB(x 



dx 



(16) 



\h [xuh{vo)]\ ~ 7= * pe < = \k [x uh (uj ] 
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where n(x), B(x) are plasma density and magnetic field correspondingly. The decay 
points Xi, X2 are situated in the vicinity of pump wave UHR resonance xuh{wq). It 
can be shown that for real plasma parameters the distance between UH resonances 
xuh{ui) — xjjh(cu ) is substantial to provide 

This allows us to neglect fa(xi t2 ) in the decay condition writing it as 

k (xi,2) = k 2 (xi,2) (17) 

To solve this equation we assume the dielectric permeability e to vary linearly in the 
region considered 

. X - X UH (uJo) , . X- X LH {i02) 

£{X, uJo) = -j-r j-^rr, e(x, u 2 ) = -yj -, — rr- 

L [x UH {uJo)\ L [x L h{w2)\ 

and obtain from ([T0jl . ffTTT . (fTTjl 



x 1 



fc (xi) = k 2 (x l ) = k (x 2 ) = k 2 (x 2 ) 



UH 



where following parameters are introduced 
Xlxt + X 3 nk A 



~< = ^^tF - < X l = L [zunM] XI = L [x LH (u 2 )] 4( W2 ) 

and x UH denotes the largest solution of the equation 

K x UH (u ) - x lh {uj 2 ) 



X UH A 2 — A 



Equations (fTHj) determine in the indistinct form decay point positions in question. 
They allow us to calculate the parameters necessary for formulation of threshold power 
equation. In particular, the phase, gained in the feedback loop, can be represented as 



$ = 77 (fa - fay + - ~ -A 



7T 



3 2 3 ^uh 2 

The length of the coherence region can be determined as 

g2{xi2) ^ 2L(x 1 , 2 V , 1 (x,,) (ig) 

Last important parameter is the value of LH wave group velocity in the decay 
points. It can be estimated as 



-V 2 (X!) ~ v 2 {x 2 ) ~ r T \u 2 



„2 t s 2 ^2^ 



2 V 2 



1/3 



A 2 u;2 e 

We will be interested in the absolute instability threshold for mode $ = 2ir. This 
mode has apparently almost the same threshold as fundamental mode N = 0, which 
has the lowest one, but still can be described in WKB approximation. In this case 

\z{ Xl )\ ~ \z{x 2 )\ ~ -eMifa^MviM 
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and to estimate the threshold we should solve an equation 
27re^l z l 
\Z\ \T(i\Z\)\ 2 ~ 

which gives \Z\ ~ 0.110. Substituting obtained expressions for decay points ()18|) and 
coherence region length (|19|) to (J13|) one obtains, that \Z\ should be calculated as 



L 4/3 



2 V 2 



-1/3 



8^ 
P 2 



8ttP(uJq — u^/Iuq), where P is the 



,4mJ ^ ceUo e T {u )4\u 2 ) 

To obtain the threshold we have taken here | a-o | 2 : 
pump wave power per unit square (in erg/(s • cm 2 )). 

Taking into account that for typical ST parameters in UHR uq ~ u pe , and 
considering maximum k z ~ ^/(SVre), when we can neglect LH wave Landau damping, 
we obtain for $ = 2tt a scaling for the threshold power P* 



P* [W/cm 2 ] = 1.4 ■ 10" 



W 



yl/3jil3/6^i/3 



.cm^GHz^eV^T 1 / 3 . 
where /o[GHz] is the probing frequency, T e [eV] is the electron temperature, L[cm] is the 
density inhomogeneity scale (|16|). B[T] is the magnetic field in plasma. 

We calculate P* for MAST tokamak parameters: fo = uJo/(2ir) = 57.5 GHz, 
T e = 100 eV, B = 3.2 kGs (in UHR position), L = 5 cm. In this case one obtains 
P* = 0.9 MW/m 2 , which gives for gaussian antenna beam with radius p = 10 cm 
threshold power P* ~ 30 kW. 



4. Discussion 



At first we discuss the approximations used. Our analysis is performed in WKB 
approximation, which holds true, when two following conditions are satisfied: 

• Decay points x\ t 2 are situated far enough from LH wave turning point x*. More 
accurately, taking into account that electric field of LH wave in the vicinity of the 
turning point can be expressed in terms of Airy function, it can be written as 

xi, 2 - x* > i A (20) 
where Airy scale £a = 2 2 ^ 3 \ 2 . In our case 

2/3 



Xl.2 



4 V 2 

and the condition (j2*U|) can be shown to be satisfied even for $ = 2ir. 

Length of decay region is not larger than Airy scale, which provides that all decay 

region is situated far enough from the turning point. The coherence region size (JT§j) 



can be estimated as 



(2AQX* 



.1/2 



and it can be shown that condition 



can be satisfied for wide range ST experiment parameters. 



(21) 
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Figure 2. Approximate representation of the dispersion curves in the vicinity of 
decay points. 



• The distance between extraordinary wave cut-off and UHR, which can be estimated 

as 

a Uce T 
Ax ~ L 

should be much larger than pump wavelength in the decay region. The last can be 
estimated as Ao(xi,2) — 2ir/x*. Corresponding condition, which provides WKB- 
representation of the UH waves to be correct, takes the form 
= u ce (x UH )Lx* i 
^ 27rio>o 

This criterion is rather strict due to low magnetic field, which is typical for ST. But 
it can be shown to be satisfied for MAST experiment parameters, where \i ~ 6. 

Our consideration, which is based on the formulae (|14|). seems to be sensitive to 
the possibility to consider decay points as separate amplifiers of the incident wave. The 
condition for that is x 2 — x\^> £, which is equivalent to 



x 2 - Xi 



2 V 2 



lV3 p 

> 1 (22) 



A 2 

Comparing that with (j2lj) . one obtains that (}2*2*j) can be satisfied for $ = 2n in rather 
narrow range of parameters. But, actually, an accurate analysis shows, that for the 
dispersion curves behavior, which in the region in question can be approximated as 
in figure |2l the decay points joint influence is the same as given by our consideration. 



5. Conclusion 



In the paper absolute instability of parametrical reflection of upper hybrid wave is 
analyzed in WKB approximation. The reflection is assumed to be accompanied by 
radiation of lower hybrid wave. Equations, describing the decay, obtained in kp- 
approximation. The decay threshold is determined. It is shown that threshold power 
increases with pump frequency and electron temperature. Threshold power is estimated 
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for typical parameters of experiment in MAST tokamak. It is shown that in this case 
parametrical reflection arises, if probing power exceeds 90 W/cm 2 , which gives 30 kW 
in a beam of 10 cm radius. 
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